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Abstract 

An optimization-based integrated controls-structures design methodology for a class 
of flexible space structures is described, and the phase-0 Controls-Structures-Integration 
evolutionary model, a laboratory testbed at NASA Langley, is redesigned using this 
integrated design methodology. The integrated controls-structures design is posed as a 
nonlinear programming problem to minimize the control effort required to maintain a 
specified line-of-sight pointing performance, under persistent white noise disturbance. Static 
and dynamic dissipative control strategies are employed for feedback control, and parameters 
of these controllers are considered as the control design variables. Sizes of strut elements 
in various sections of the GEM are used as the structural design variables. Design guides 
for the struts are developed and employed in the integrated design process, to ensure that 
the redesigned structure can be effectively fabricated. The superiority of the integrated 
design methodology over the conventional design approach is demonstrated analytically by 
observing a significant reduction in the average control power needed to maintain specified 
pointing performance with the integrated design approach. 
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I. Introduction 


Currently spacecraft are designed in separate iterative sequences within the structural 
and control disciplines. The structural design, performed first, takes into account loading 
considerations that would occur during launch, reboost, or operational maneuvers. The sizes 
and masses of mission-related components are estimated and a structure is designed which 
maintains the desired spatial relationships among various components during operations. 
In the next stage, a controller is designed for the fixed structure to orient, guide, and 
move the spacecraft to obtain the required spacecraft performance. Measures of spacecraft 
performance may take different forms such as limits on pointing jitter, transient response 
deviations, control power requirements, and so on. The control design must also provide 
satisfactory closed-loop stability robustness. This separate-discipline approach has been 
used successfully in most of the missions in the past. It works well when a relatively high- 
stiffness structure is feasible, the nonstructural components are concentrated masses and 
inertias; and the performance requirements of the missions are not stringent. In these cases 
the structural modes are beyond the controlled bandwidth, so that there is minimal control- 
structure interaction. However, this approach will not meet, the stringent performance 
requirements of future space structures. A number of future space missions will utilize 
flexible structures in low-Earth orbits and geostationary orbits. Examples of such missions 
include space science platforms, space processing facilities, and Earth observation systems. 
Such missions typically require distributed-ma.ss components such as booms, solar arrays, 
and antennas, whose dimensions range from a few meters to possibly hundreds of meters. 
To minimize the costs of construction, launching, and operations, it is necessary to make 
the structure as light as possible. However, the combination of large size and low structural 
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mass leads to increased flexibility and makes attitude control of the structure to a specified 
precision more difficult. 

12 

Controls-structures interaction (CSI) in the form of destabilizing spillover has been 
verified in simple Earth-based laboratory experiments as well as in the design, analysis, 
ground development, test, and flight operation of space systems in industry^. The current 
approach to solving CSI problems is to design the spacecraft to avoid undesired dynamical 
interaction. This generally requires either stiffening the structure or slowing down the 
control system response. Stiffening the structure simplifies the control design problem, in 
that the predominant dynamics tend toward rigid body, but is costly in terms of mass as 
well as launch packaging, leading to increased fuel consumption. Slowing down the control 
response produces control inputs which have less of a chance of producing destabilizing 
effects, but is costly in terms of reduced performance capability. Neither approach is 
completely satisfactory. What is needed is a new design approach which has the capability 
of avoiding any damaging aspects of controls-structures interaction, while at the same time 
identifying and exploiting the beneficial aspects. 

A number of studies have been presented in the past, decade using integrated control- 
structure design methodologies to enhance the overall performance of the controlled system. 
Minimization of structural mass, subject to constraints on closed-loop eigenvalues, damping 
parameters, and structural natural frequencies, has been considered in Refs. 3 and 4. While 
dependence: of control cost, on the structural parameters is clear, since system matrices 
are functions of structural parameters, total structural mass is not. considered directly 
as a. function of control parameters, or control effort., in these studies. The effect, of 
variable actuator masses, as a. function of control effort., is addressed in Refs. 5 and 6. 
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Other studies ’ have considered minimization of a control objective such as the norm 
of feedback gain matrices, as a measure of control effort, with constraints on closed-loop 
responses and total structural mass. A covariance approach is used in Ref. 8 to constrain 
overall system response with bounds on the allowable control energy. Multiobjective 
optimization has been used to optimize a combination of structural objectives such as 
structural mass, and control objectives, such as control energy or quadratic objective 
functions obtained from linear regulator theory 9-11 . Multiobjective optimization provides 
tradeoffs between competing objectives, thus providing a continuum of system designs, with 
different performance characteristics. Reference 12 includes a good survey of numerous 
other approaches employed for integrated control-structure design. However, most of the 
techniques in the literature have been applied to simple analytical models or laboratory 
apparatuses. Only recently have there been some studies in integrated design of large-order 
flexible space structures (with thousands of degrees of freedom) 5,5,15,11 . 

The OSI program at NASA Langley 15 has been developing an integrated controls- 
structures design methodology which is described in this paper. The design methodology 
is based on the recognition that a high degree of coupling exists between the control and 
structural disciplines in the control of flexible space structures. Rather than performing 
Separate structural and control designs, a unified environment for integrated control- 
structural modelling, analysis and synthesis is developed. Within this environment, a design 
iteration consists of updating all critical (control and structure) design variables in a single 
integrated computational framework by assessing controlled performance with the current 
design variables. Parameter optimization techniques, employing nonlinear mathematical 
programming, are used for the synthesis of an optimal integrated structure. Realistic 
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constraints, which have not been considered in previous literature, were imposed on the 


current design process, to ensure that the structure obtained through integrated design 
could be practically fabricated. For example, the strut designs were obtained using design 
curves relating effective areas and effective densities that imposed the requirements of 
manufacturability on the design process. This structure was subsequently used to verify 
the benefits of control/structure integrated design experimentally, as described in Ref. 16. 

The phase-0 CSI evolutionary model (CEM), a laboratory test structure at NASA 
Langley (shown in Fig. 1), was used as the candidate flexible structure for controls- 
structures integrated redesign. Two dissipative control strategies, namely, static dissipative 
and dynamic dissipative controllers, were employed for feedback compensation, since these 
compensators guarantee closed-loop stability in the presence of unmodelled dynamics and 

. . . 17 

parametric uncertainties . The performance measure for optimization is the average control 
energy required to maintain prescribed line-of-sight. pointing, in response to persistent, 
white noise at specified disturbance locations. Integrated redesign of the phase-0 CEM was 
performed with the two control strategies to improve the performance of th^r'- controlled 
system. Concurrently to the integrated design process, control-optimized dissipative 
controllers were designed for the nominal phase-0 CEM. A comparison of the controlled 
performance of the integrated designs with that of the conventional control-optimized 
designs indicate over 40% reduction in control power, while providing the same pointing 
performance. These results clearly demonstrate the advantage of integrated controls- 
structures design methodology over the traditional sequential approach. The redesigned 
structure was fabricated and assembled in the laboratory and a number of active control 
experiments were performed with it to experimentally validate these results. This effort 
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represents the first experimental validation of the integrated design methodology for realistic 
flexible space structures, and is described in Ref. 16. 

II. Mathematical Models 

Structural Model 

The dynamics of flexible space structures can be treated as dynamics of a distributed 
parameter system from a continuum mechanics viewpoint. Applying the principles of 
classical mechanics to the distributed parameter system results in a system of partial 
differential equations with associated initial and boundary conditions. However, this 
framework for description of structural dynamics is of limited value in characterizing the 
dynamics of realistic, flexible-Space structures. Space structures are conventionally analyzed 
using approximate 1 tffjhniques such as the finite element method. A finite-dimensional, 
linear, time-invariant, mathematical model of a flexible space structure is given as 

Mx + Dx + Kx = B ii + Hw (1) 

where X is an n X 1 displacement vector; U is the m X 1 control input vector; w is the 
p X 1 disturbance vector; M is the positive-definite inertia matrix; D is the positive serni- 
definite open-loop damping matrix; A is the positive semi-definite stiffness matrix; B is an 
n X m control influence matrix; and H is n X p disturbance influence matrix. The rate 
and performance output equations are 

Ur = (-'rX ! Uper = ( 'per X (2) 

where y r is the rate measurement; y 1}er is a performance vector; and C r , and C p er are 
the rate output and performance output influence matrices, respectively. The second-order 
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representation of the structure given in Eq. (1) is obtained by finite element modeling of the 
structure, and the order of these models can be quite large. For design and analysis purposes, 
the order of the system is reduced using a modal truncation approach, wherein significant 
modes in the input-output, characterization of the plant are retained and the remaining 
modes are truncated. The system equations in modal coordinates for the retained modes 
are written as 

.1 / , (j , T DfCjf T /\ , (j , — $ Bu T $ Hw 

A) 

= I r u T T d lu 

where q r is an r X 1 vector of modal amplitudes; M r , I ) , . A r are, respectively, the 
generalized inertia, damping and stiffness matrices; and $ is an n X r matrix whose columns 
are the r structural eigenvectors associated with the included modes. If the mode shapes 
are normalized with respect to the inertia matrix, and modal damping is assumed, then 
M r = I r x r i D r = Diag{2(iiOi , . . . , 2( r iO r }, and K r = Diag{ujf, . . . , where 
LOi and (,j, l = 1, . . . , r are the open-loop frequencies and damping ratios. With collocated 
and compatible measurement sensors and control actuators, 

y r = B T $q r = T c q r (4) 

Note that the collocation of these sensors and actuators is necessary for the implementation 
of the dissipative controllers. This collocation guarantees the system to be minimum phase, 
i.e, transmission zeros and poles of the system are in the left-half plane, and enhances the 

stability robustness of the overall system" . 

f 7 1 * 7 1 1 J 1 

Defining the state vector, z = \q^ qg: } , the dynamics of the system can be written 
in a first-order form, 

z = Az + Bu + Hw (5) 
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where 


A = 


0 Irxr 

— Kr —Dr 


; B = 


r T 

c 


; H = 


rT 

L 1 d J 


The rate output and performance output vectors are given by 


y : — r c . — ( r . . Uper — ^'per^ ^ pi i - 


Here, C per is the output influence matrix associated with the performance output, y pe 


Controller Design Methods 

Control system design for flexible space structures is a challenging problem because of 
their special dynamic characteristics, which include a large number of structural modes 
within the controller bandwidth; low, closely-spaced structural frequencies; very small 
inherent damping; and a lack of accurate knowledge of modal parameters. Moreover, 
the controller must be of a reasonably low order to be implementable; and must also 
satisfy the performance specifications (i.e., constraints on root.-nrean-square (RMS) pointing 
error, or the desired closed-loop bandwidth). The controller must also have robustness 
to nonparametric uncertainties, such as unmodeled structural modes, and to parametric 
uncertainties, that is, errors in the knowledge of the design model. 

Two major categories of controller design methods for flexible space structures are 
model-based controllers and dissipative controllers. A model-based controller generally 
consists of a state estimator (a Kalman-Bucy filter, or some other observer) followed by a 
state feedback controller. The state estimator utilizes the knowledge of the design model in 
its prediction part. Using multivariable frequency-domain design methods, such controllers 


can be made robust to unmodeled structural dynamics, for example, the spillover effect 
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can be overcome". However, model-based controllers generally tend to be very sensitive 

to uncertainties in the design model, in particular, to uncertainty in the structural mode 
2 17 

frequencies"’ . An analytical explanation of this instability mechanism may be found in 
Ref. 17. 

Dissipative controllers utilize special passivity properties of structural dynamics, with 

collocated and compatible actuators and sensors, and provide robust stability in the presence 
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of both nonparametric and parametric uncertainties"’ . Therefore, dissipative controllers 
have been used for feedback control in the integrated design methodology. 

A. Static Dissipative Controller 

The simplest dissipative controller is the static or constant-gain dissipative controller. 
Using collocated and compatible actuators and sensors, such as torquers with attitude and 
attitude rate sensors, and thrusters with linear position and velocity sensors, the constant- 
gain dissipative control law is given by: 

u = G r y r G pi/p ( 6 ) 

where y r and y p are m X 1 rate and position measurement vectors, where m is the number 
of sensors ; G r and Gp are m X m symmetric, positive semi-definite rate and position gain 
matrices, respectively. This control law has been proven to give guaranteed closed-loop 
stability despite unmodeled elastic modes, parameter errors, certain types of actuator and 
sensor nonlinearities (such as saturation and dead zone), as long as actuator dynamics are 
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limited to zero-order or first-order forms". For space structures with zero frequency, rigid- 
body modes, position feedback is essential to ensure stability of the closed-loop system. 
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However, for the ground-based structure considered in this paper, with non-zero frequency 
suspension modes, rate feedback alone is sufficient to guarantee stability. Therefore, for the 
integrated design and experimental studies, only rate feedback is considered. The drawback 
of this controller is that the achievable performance is inherently limited because of its 
simple mathematical structure. 

B. Dynamic Dissipative Controller 

To to obtain higher performance while still retaining the highly desirable robust 
stability, dynamic dissipative compensators can be used. An ?? c -order, t.wo-level, dynamic 
dissipative controller is given by: 


x c — A c x c T B c y r ( i ) 

u = C c x c G r y r Ct py.p ( 8 ) 

where A c , ll< . and C c are the compensator system, input, and output matrices, respectively. 
G r and Gp are symmetric, positive semi-definite rate and position gain matrices correspond- 
ing to the static dissipative inner loop. For ground test articles, with no zero frequency rigid- 
body modes, position feedback is not necessary. The dynamic compensator, in the outer 
loop, is dissipative if its transfer function is strictly positive real. The positive realness 
lemma, or the Kalman-Yacubovich lemma , ^ assures this condition when A c is Hurwit.z, 

(A C7 B C ) is controllable, (A c ,C c ) is observable and there exist a matrix P = P ^ > 0 
such that 


A t c P + PA c = -Q 


c c = b t c p 


( 9 ) 

( 10 ) 
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is satisfied, where Q = Q ^ > 0. The main characteristic of dissipative controllers is 
that they do not rely on the knowledge of the design model to ensure stability. However, 
these controllers do use this knowledge to obtain the best possible performance. This 

controller assures robust asymptotic stability regardless of unmodeled structural dynamics 

1 7 

or parametric uncertainties 

III. INTEGRATED DESIGN OF THE CEM 

This section presents the integrated design methodology, and its application to the 
redesign of the nominal phase-0 OEM. 

The phase-0 OEM, shown in Fig. 1, is a laboratory testbed at NASA Langley, for 
experiments in control of flexible space structures. It consists of a 62-bay main truss (each 
bay is 10 inches long), two vertical towers, and two horizontal booms. The structure is 
suspended from the ceiling about. 840 inches above the main truss using two cables attached 
to the horizontal booms, as shown in Fig. 1. A laser source is mounted at the top of the taller 
tower, referred to as the laser tower, and a reflector with a mirrored surface is mounted on 
the other tower, called the reflector tower. A laser beam is reflected by the mirrored surface 
onto a detector surface 660 inches above the reflector. The line-of-sight. pointing problem is 
to maintain the laser close to its nominal position on the detector surface, despite persistent 
disturbances. Eight proportional bidirectional air thrusters, with maximum output force of 
4.4 lbs each, are available at stations 1 to 8, shown in Fig. 1. Almost collocated with the 
thrusters at the eight stations are servo acjteleromet.ers to provide output measurements. 
Accelerometer signals are subsequently integrated with the aid of wash-out filters to provide 
the required rate information. A detailed description of the pha.se-0 CEM is available in 
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Ref. 19. 


The phase-0 CEM is representative of flexible space structures in that a number of low 
frequency closely spaced modes are within the bandwidth of its controllers, inherent damping 
is low, and its modal parameters are uncertain. However, one major difference between any 
ground-based structure, which must be suspended under the influence of gravity, and space 
structures is the absence of zero frequency rigid body modes in ground-based structures. 
The finite element model of the system has .3216 degrees of freedom; therefore, a large 
computational effort is required for solution of the structural eigenvalue problem of that sizq. 
The control design model consisted of the first 30 modes of the structure, which includes 
24 flexible modes and six suspension modes (non-zero frequency, rigid-body modes, due to 
suspension of the structure in the presence of gravity). A modal damping ratio of 0.1% 
was assumed. The modal frequencies (in Hz) of the first 10 modes of the nominal phase-0 

CEM are {0.147, 0.149, 0.155, 0.730, 0.748, 0.874, 1.473, 1.738, 1.882, 2.294}. The 

first six modes, ranging from 0.147 Hz to 0.874 Hz, are the suspension modes. Modes 7 
and 8 are the first two lateral and vertical bending modes, and mode 9 is the first torsional 
mode of this structure. 

For the integrated design problem, white-noise disturbances of unit intensity are applied 
to the structure at stations 1 and 2, and the feedback control inputs are applied at 
stations 3 through 8. Past experience with the phase-0 CEM structure had shown that 
actuators at stations 7 and 8 easily destabilize the laser tower modes, in the form of 
spillover destabilization, which is typical in control of flexible structures. Therefore,: even 
though actuators 7 and 8 are very efficient locations for exciting the structure, they were 
included for feedback control in order to consider the robust stability issue objectively. After 
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actuators 7 and 8, actuators 1 and 2 are the most effective actuators in exciting the structure. 
Therefore, actuators 1 and 2 were chosen as disturbance sources for integrated design and 
experimental validation. For the design process, no sensor and actuator dynamics were 
included. Thus, the system dynamics are given by Eq. (5), where the matrix r c contains 
the modal displacements at the control stations 3 through 8; and matrix contains modal 
displacements at the disturbance stations 1 and 2. The performance vector, y per , is the line- 
of-sight pointing error or deviation of the laser point on the detector system from its nominal 
position in local x and y coordinates. The performance, output matrix, C per , is computed 
from the kinematics of small angular motion about, the nominal configuration of the system. 

The feedback control configuration used for the synthesis of the integrated cont.rols- 
structures design as well as the control-optimized design is as follows. Persistent disturbance 
noise, w(t), is applied to the structure at the disturbance locations (stations 1 and 2), and 
the deviation of the line-of-sight. pointing error vector, y per (t), as measured by the laser- 
detector, is to be maintained within desired specifications. The feedback control system 
applies control inputs, u(t), at control stations 3-8, using the measured outputs, y r (t), at 
nearly collocated nodes. 

The design optimization problem is to minimize the steady-state average control power 
while maintaining a specified root.-mean-square (RMS) line-of-sight pointing error, and 
without exceeding the total mass budget of the nominal pha.se-0 CEM. Mathematically, 
the design problem is to minimize 


J = lim \Tr(E{u(t)'U T (t)})] 

t — .-OG 


(ID 
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with respect to structural and control design variables, subject to the constraints: 

lim [Tr{E{y per {t)y^ er {t)})] 2 < M str < M max (12) 

t — '-OO 

where 1 is the maximum allowable line-of-sight pointing error, M max is the mass 
budget, Tr denotes trace of the matrix, and E is the expectation operator. 

For the redesign of the CEM, the limit on the root.-mean-square line-of-sight pointing 
was chosen as y'”"/ equal to 2.4 in, which is an order of magnitude reduction from the open- 
loop RMS line-of-sight pointing for phase-0 CEM structure (which is 22.54 in); and M max 
was chosen to be 1.92 lb. sec" /in, which is the nominal mass of phase-0 CEM structure. 

Using the static or dynamic dissipative controller, the closed-loop system dynamics can 
be written as 

x = Ar + II ir\ y per = C per x; u = C u x (1.4) 

where X is the state vector for the closed-loop dynamics, w(t) is zero mean, white noise 

disturbance applied at the disturbance locations, y.p er is the line-of-sight pointing error, u(t) 

is the control vector, and A, / / . C per and C u are the corresponding closed-loop system 

matrices. The steady state covariance matrix for the closed loop state, is computed by 

9D 

solving the following Lyapunov equation 

AEj + Ti x A T + HTi w H T = 0 (14) 

where is the covariance matrix for the disturbance noise, w(t). The steady-state average 
control power is given as 

lim Tr{E{u(t)u T (t)}} = Tr[C. u T, x C^] (15) 

t — ^OG 
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and root-mean-square line-of-sight pointing error is 


lim [Tr(E{y per (t)y^ er (t )})] 2 = [Tr(C pe (16) 

The total mass of the structure, M s f r , is obtained from the mass matrix of the finite 
element analysis routines. 

For the structural design, the CEM structure was divided into seven sections, shown 
in Fig. 2, namely, three sections in the main truss, one section for the laser tower, one 
section for the reflector tower, and one section for the two horizontal booms. The primary 
dynamic characteristics of these sections were defined by the mass and stiffness properties 
of the struts and node balls comprising these sections. A typical strut and node ball of the 
resigned CEM structure is shown in Fig. .3. The strut/node ball system was modeled as an 
axial rod between node ball centers representing the stiffness of the strut and node ball with 
a point mass at each end representing the mass of the node ball. The stiffness of the axial 
rods in each section were used as the structural design variables. From a modeling point 
of view the rod stiffness is defined by its effective area, length, and modulus. Civen that 
the length and modulus of the struts were held constant, the effective area of the axial rod 
was used to quantify the stiffness of the struts. Therefore, three structural design variables 
were used in each section, namely, the effective areas of the longerons, battfns and diagonals. 
Thus, the integrated design of the structure involved a total of 21 structural design variables. 

In an ideal design with uniform struts, the effective area of the tod would be the cross- 
sectional area of the strut. However, in this case the strut is a complex member with joints 
and variable cross-sections. What is required is a means to relate the effective area used to 
model the strut/node ball system to a strut that can be manufactured. The strut of Fig. 3 
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is composed of three sections, namely, the node ball, tube, and connection hardware. The 


system is designed such that the node ball and connection hardware remain the same for 
all struts. However the tube portion of the’ strut can be changed. From a stiffness point 
of view the strut. /node ball system can be thought of as three springs in series. The end 
springs represent the stiffness of half a node ball and connection hardware, while the center 
spring is the stiffness of the tube. Thdeeffective area can be represented as: 

_ {k c k t )£ 

' ff (2 h + k c )E 

where k c is the stiffness of the end springs, kf is the stiffness of the tube, £ is the distance 
between node ball centers, and E is the modulus of elasticity. If the center -to-center length, 
modulus, and stiffness of the node ball and connection hardware remain constant, then the 
effective area will be governed by the stiffness of the center tube. Associated with each 
effective area is an effective density defined as 


Peff 


A eff? 


where m is the total mass of the strut, excluding the node ball mass. The effective density 
combined with the node ball mass is used to define the mass of the strut./ node ball system 
in the finite element, analysis of the structure. 

Now the stiffness and mass properties of the structure are functions of the geometry 
of the center tube of the strut.. To ensure that, the strut. /node ball system could be 
manufactured, design guides were semi-empirica.lly developed that, defined a. two dimensional 
design space of manufacturable struts" . The design space gives the relationship between 
the effective area, and effective density for the strut, design shown in Fig. .3. The strut, design 
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guides are shown in Fig. 4, where the dotted line corresponds to longerons and battons, 
and the solid line is the design guide for diagonals. The curve represents the minimum mass 
(density) strut for a given strut stiffness (effective area). The curve is continuous due to 
the manufacturing process. The center tube is manufactured by turning a stock tube to 
a specified outer diameter, thereby, allowing a variety of center tube stiffness. The design 
space is limited at low effective areas (towards the left end of the design guides in Fig. 4) 
by the strength of the strut, whereas the high effective area (towards the right end in Fig. 
4) is limited by the linearity of the connection hardware. The control design variables for 
static and dynamic dissipative controllers are described in the next section. 

An integrated design software tool, called GSI-DESIGN, being developed at the NASA 
Langley, was used to perform the parameter optimizations numerically. G.SI-DESIGN 

OO nA 

uses an in-core database architecture and is composed of public domain software"" " . 
The package has control, structural and optimization modules linked together in a unified 
environment to perform design iterations on both structural and control design variables. 
A description of the contents of the GSI-DESIGN structural module may be found in Ref. 
25. Employing a four-processor Alliant FX-80 digital computer, integrated optimization 
was performed using the Automated Design Synthesis (ADS) software"". This design 
optimization took about. 8 hours of cpu time for each run. Gradient computations were 
performed using finite difference approximations. An interior penalty function method 
of ADS was used to solve the nonlinear programming problems. In this method, the 
constrained optimization problem is transformed into an unconstrained problem through 
creation of a pseudo-objective function, which is the sum of the original objective function 
and an imposed penalty function (a function of the constraints"®). The Reverse-Cut.hill- 
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27 

McKee algorithm" for minimizing the bandwidth of the banded stiffness and mass matrices 
was used to reduce computational and memory requirements. Additionally, analytical 
expressions for eigenvalue/eigenvector sensitivity, with respect to the structural design 
variables" , were used in the integrated design process to approximate the eigenvalues and 
eigenvectors at design points which are in the neighborhood of the nominal design point. 
This approximation was in the form of a first-order Taylor series approximation and resulted 
in substantial computational savings since it. removes the need for costly computation of 
structural eigenvalues and eigenvectors at many of the optimization moves. The Lyapunov 
equation for closed-loop state covariance matrix, Eq. (14), is solved by transforming the 
matrices to a Schur form and solving the resulting linear equations, as described in Refs. 2.3 
and 29. The control design variables used for static and dynamic dissipative control, and 
the results of the integrated design process are described next. 

IV. Integrated Design Results 
A. Integrated Design with Static Dissipative Controller 

This section describes the control design variables for static dissipative controllers and 
the design optimization results. Since the test article does not have any zero frequency rigid 
body modes, only rate feedback is employed for the static dissipative controller, i.e. 

u = Cr r y r , (17) 

where the rate gain matrix, G r , is chosen to be a 6 X 6 diagonal matrix, whose; elements 
constitute six control design variables. Note that a diagonal gain matrix is used for 
simplicity, and because the resulting decentralized controller generally exhibits superior 
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performance robustness. The closed-loop matrices in Eq. (13) are 


A 


0 1 rxr 

—K r —D r — TjG r T c 



— [o <s'. r r c ]. (is) 


Average control power and root.-mean-square line-of-sight pointing error for static dissipative 
controllers, are computed using Eqs. ( 1 5 ) - ( 16). A total of 27 design variables, 21 
structural design variables and 6 control design variables, were used for the integrated 
design optimization for the static dissipative integrated design. 

The results of the design optimizations are summarized in Table 1. First, the average 
control power was minimized with respect to the control design variables, while the structural 
design variables are fixed at the nominal phase-0 values, namely, effective areas for battons 
and longerons being 0.134 in" and effective area for diagonals being 0.124 in". This 
corresponds to the conventional approach of synthesizing the best achievable controllers 
for a fixed structure. The control-optimized design required an average steady-state control 
power of 7.11 lb" to maintain RMS line-of-sight error within the specified limit, that is, 2.4 
in. Next, an integrated design was performed to minimize the average control powqr- with 
respect to both control and structural design variables. The results (see Table 1) indicate 
an average control power of 4.21 lb" to maintain the same RMS line-of-sight error. The 
integrated design results in a reduction of more than 40% in the average control power over 
the conventional design, for the same RMS line-of-sight pointing performance. The effective 
areas for the structural design are shown in Table 2. Comparing the new design variables 
with the effective areas of the nominal phase-0 CEM, it is observed that the longerons 
for all three sections of the main truss, particularly the section closest to the disturbance 
sources, are considerably stiffened. The longerons in the laser tower are also stiffened. On 
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the contrary, the horizontal booms and the reflector tower became more flexible, partially to 


satisfy the mass constraint. Generally, all the diagonals and the battens decreased in size, to 
satisfy a constraint on the total mass i.e. the mass of the redesigned structure had to be less 
than or equal to the mass of the nominal phase-0 design. Consequently, mass was taken from 
the battens and diagonals and was redistributed to the longerons of some sections because 
longerons are quite effective in increasing the stiffness of that section. This trend may be 
attributed to a trade off between structural controllability and observability at the actuator 
stations, structural observability of the performance vector and structural excitability at the 
disturbance locations. The areas near the disturbance sources, i.e. stations 1 and 2, were 
stiffened in order to reduce the sensitivity of the structure to external disturbances at those 
locations, while ensuring tha|; no appreciable loss of controllability and/or observability 
occurred at the control stations. The control gains for the control-optimized and the 
integrated designs are shown in Table 3. Generally it is observed that the gains for the 
redesigned system are considerably lower than those for phase-0 GEM design (except for 
actuator 6). This is to be expected since the required control power for redesigned system 
is significantly lower than that of the nominal phase-0 GEM. 


B. Integrated Design with Dynamic Dissipative Controller 

The dynamic dissipative controller represented by Eqns. (7) and (8), with no static 
inner loop, that is, G p = 0 and G r = 0 , was used for the integrated redesign of the 
structure. The compensator state matrix A c and compensator input influence matrix B c , 
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were assumed to be block diagonal, as shown below, where Ap and Bp ■ l = 1, 2, . . . , 6 
are, respectively, 2x2 matrices and 2x1 vectors, defined as 

n li roi 

Ap = n i B c i = .. (21) 

Qj pi 1 

The parameters, Ct i and f3{, for i = 1,2, . . . , 6 , are the control design variables. Further- 
more, the weighting matrix Q in Eq. (9) was assumed to be diagonal, that is, 

Q = Di ag ( <?l , t/ 2 , • • • , t/ 12 ) , (22) 

with Cj j , for j = 1,2, ..., 12, also being the control design variables. For the dynamic 
dissipative controller, the closed-loop system matrices in Eq. (13) are 

_ r A —BC c 1 \H] _ r 

- 4 =Ug .4, J ; r '=[ r '- °] < 23 > 

For average control power computation, C u = — [0 ■ Again, Eqns. ( 15 )-( 16 ) are used 

for computing control power and root-mean-square line-of-sight pointing error. Thus, the 
scalar variables a i, ftipi = 1, 2, . . . , 6 and C/j,j = 1, 2, . . . , 12 constitute 24 control 
design variables. Alongwith 21 structural design variables, the total number of design 
variables for integrated redesign with dynamic dissipative controllers is forty five. 

Table 1 shows the results of designs with the dynamic dissipative controller. The 
control-optimized design for phase-0 CEM required a control power of 6.41 lb “ to maintain 
RMS line-of-sight pointing within 2.4 in. The integrated redesign reduces the average 
control power by 44% over the conventional design to 3.64 lb 2 . The effective areas for 
the structural elements for integrated design with the dynamic dissipative compensator are 
shown in Table 2. Effective areas for longerons, battens and diagonals for the integrated 
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design with dynamic dissipative controllers show the same trends as that for the static 
dissipative controller. Control design variables for the control-optimized design variables 
and the integrated design are shown in Table 4. Figure 5 compares the maximum and 
minimum singular value plots of optimal dynamic dissipative controllers for phase-0 CEM (in 
dashed lines) and redesigned system (in solid lines). It is observed that the controller gains 
are generally smaller throughout the frequency spectrum except in the very low frequency 
region. 

The results obtained for both the static and dynamic dissipative controllers clearly show 
that integrated controls-structures design methodology can yield a substantially superior 
overall design over the conventional sequential design scenario. Moreover, Table 2 shows 
that the optimal structure for both dissipative control laws exhibit similar trends. In fact, 
the structural design variables for the two optimal structures are within 20% of each other. 
A structural design close to both these designs was chosen for fabrication and assembly in the 
laboratory. This redesigned structure was referred to as phase- 1 CEM. A number of active 
control experiments were performed with both the structures, that is, phase-0 CEM and 
phase- 1 CEM, to validate the benefits of integrated controls-structures design demonstrated 
analytically in this paper. Thistftxperimental validation is described in Ref. 16. 

V. Concluding Remarks 

An optimization-based integrated controls-structures design methodology has been 
presented in this paper, using static and dynamic dissipative compensators. To demonstrate 
the benefits of this design approach, the nominal phase-0 CEM structure was redesigned, 
using the integrated design methodology, to minimize the average control power required 


22 



to maintain specified root-mean-square line-of-sight pointing error under persistent distur- 
bances. The redesign with static dissipative compensators resulted in a 40% reduction in 
control power, whereas the redesign with dynamic dissipative compensators produced a 
44%) reduction in control power. Therefore, it has been demonstrated, analytically, that 
integrated controls-structures design can yield designs which arfe.substantially superior to 
those obtained through the traditional sequential approach. Moreover, in order to show that 
the benefits achieved analytically can be realized in practice,, the redesigned structure was 
fabricated and assembled in the laboratory. Experimental validation of these results has 
been performed successfully, and will be presented in a subsequent paper. Finally, this work 
demonstrates the capability of a software design tool (CSI-DESIGN) which implements the 
automated design procedure in a unified environment for structural and control designs. 
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Table 1 Results of design study with static and dynamic dissipative controllers. 



RMS Line-of-sight 
Error (in) 

Control Power 
(lb"2) 

Control-Optimized Design with 
Static Dissipative Controller 

2.4 

7.11 

Integrated Design with 
Static Dissipative Controller 

2.4 

4.21 

Control-Optimized Design with 
Dynamic Dissipative Controller 

2.4 

6.41 

Integrated Design with 
Dynamic Dissipative Controller 

2.4 

3.64 
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Table 2 Structural design variables for integrated design 
with static and dynamic dissipative controllers. 



Longerons, in~2 


Battens, in~2 


Diagonals, in"2 



Var. 

No. 

Static 

Dynamic 

Var. 

No. 

Static 

Dynamic 

Var. 

No. 

Static 

Dynamic 

Sec. 

No. 

Description 

1 

0.330 

0.330 

2 

0.082 

0.077 

3 

0.082 

0.098 

1 

main truss 

4 

0.085 

0.080 

5 

0.083 

0.087 

6 

0.085 

0.087 

2 

horizontal 

boom 

7 

0.173 

0.142 

8 

0.082 

0.086 

9 

0.082 

0.082 

3 

laser tower 

10 

0.260 

0.295 

11 

0.082 

0.080 

12 

0.081 

0.066 

4 

main truss 

13 

0.257 

0.258 

14 

0.081 

0.078 

15 

0.079 

0.066 

5 

main truss 

16 

0.095 

0.100 

17 

0.081 

0.077 

18 

0.079 

0.066 

6 

horizontal 

boom 

19 

0.096 

0.117 

20 

0.081 

0.083 

21 

0.082 

0.083 

7 

reflector 

tower 
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Table 3 Control design variables for integrated 
design, with static dissipative controllers. 


Var. 

No. 

Control-Optimized 

Design 

Integrated 

Design 

Actuator 

No. 

1 

1.2110 

0.5586 

3 

2 

2.0634 

0.6837 

4 

3 

1.2985 

0.7899 

5 

4 

0.5352 

1.1117 

6 

5 

1.5050 

1.0201 

7 

6 

0.4593 

0.5982 

8 
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Table 4 Control design variables for integrated 
design, with dynamic dissipative controllers. 


Var. No. 

Control- 

Optimized 

Design 

Integrated 

Design 

Var. No. 

Control- 

Optimized 

Design 

Integrated 

Design 

A c[ 2,1] 

136.57 

135.39 


15287.63 

18122.76 

^c[2,2] 

105.25 

114.73 

Q[2,2\ 

30050.42 

12173.74 

^c[4,3] 

136.60 

301.81 

Q[3,3] 

15125.39 

10430.10 

^c[4,4] 

104.33 

102.17 

Q[ 4,4] 

29103.43 

13600.71 

^c[6,5] 

103.80 

80.42 

Q[5,5] 

19640.50 

24674.52 

^c[6,6] 

107.72 

125.07 

Q[6,6] 

24331.62 

19487.68 

C5 

r oo 

94.56 

58.21 

Q[ 7,7] 

20582.68 

22964.97 

^■c[8,8] 

121.26 

103.52 

Q[8,8] 

22159.79 

21224.58 

^■c[10,9] 

136.92 

95.28 

<9 [9, 9] 

15148.89 

22693.96 

A:[10,10] 

103.87 

107.26 

<9[10,10] 

29362.52 

20381.06 

^■c[12,ll] 

100.11 

136.11 

Q[ll,ll] 

20009.49 

17783.45 

^c[12,12] 

126.11 

133.30 

Q[12,12] 

21511.89 

15376.09 
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Fig. 1 Schematic of the CEM testbed. 
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Fig. 2 CEM displaying seven sections for integrated redesign. 
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Fig. 3 Strut design for phase- 1 CEM. 
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Fig. 4 Strut design curves, longerons and battens 
(dashed line) and diagonals (solid line). 
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Fig. 5 Maximum and minimum singular value plots 
controllers for phase-0 CEM (dashed line) and Pha: 



